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Abstract. We investigate entanglement spectra of the SO(3) bilinear-biquadratic 
spin-1 chain, a model with phases exhibiting spontaneous symmetry breaking (both 
translation and spin rotation), points of enlarged symmetry, and a symmetry-protected 
topological phase (the Haldane phase). Our analysis reveals how these hallmark 
features are manifested in the entanglement spectra, and highlights the versatility 
of entanglement spectra as a tool to study one-dimensional quantum systems via small 
finite size realizations. 


1. Introduction 

The notion of entanglement is central to quantum information theory, and has been 
applied with great success in the analysis of condensed matter systems [1]. Von Neumann 
entanglement entropy and its generalizations have been invoked in various contexts, 
including thermalization, diagnosis of topological phases, and properties of conformally 
invariant quantum systems. Beyond entanglement entropy, the full ‘entanglement 
spectrum’ of a system’s reduced density matrix p A reveals much detailed information 
about the structure of its ground state [2] . In particular, block-diagonal structures in p A 
derived from a many-body ground state projector |T 0 ) (T 0 | may be exploited to analyze 
entanglement spectra as a function of quantum numbers associated with symmetries 
which are preserved by the entanglement cut. The notion of an ‘entanglement gap’ was 
established [2, 3] to suggest a separation of “relevant” low-lying entanglement levels 
from “spurious” low probability (high entanglement energy) levels. In some cases, 
such as the fractional quantum Hall effect and symmetry-protected topological phases 
(SPTs) [3, 4, 5, 6, 7, 8, 9, 10, 11, 12], the relevant low-lying levels have similarities to 
the energy levels of an edge Hamiltonian, although this line of reasoning is not universal 
[13], 

While the task of connecting entanglement measures to physical observables is still 
at an early stage [14], the entanglement spectrum in particular has become a major 
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tool to investigate one-dimensional quantum many-body systems. On a qualitative 
level, this is due to the fact that entanglement spectra allow the reliable extraction 
of information far beyond that obtained from the entanglement entropy. For gapless 
1-d systems, the entanglement spectral distribution is fully characterized by the central 
charge of the associated conformal field theory [15] which one could deduce from the 
entanglement entropy [16]. The entanglement spectrum, however, also determines the 
boson compactification radius associated with the operator content of the boundary 
CFT [17]. For gapped 1-d systems the entanglement spectrum can completely determine 
whether a system is topological or not, and provides a classification of such topological 
1-d phases [11, 18]. 

In this article we investigate the finite size structure of the entanglement spectrum 
of the spin-1 bilinear biquadratic model for the spatial ground state bipartition of a 
12 site realization, which will be defined in Section 2. For the ferromagnetic phase 
exhibiting spontaneous spin rotation symmetry breaking in the thermodynamic limit, 
the entanglement spectrum is calculated analytically and interpreted with respect to 
finite vs. infinite size entanglement in Section 3. In Section 4, the Uimin-Lai-Sutherland 
point of the quadrupolar phase is shown to exhibit an extensive number of zeroes in 
the reduced density matrix which we analytically trace back to its enlarged internal 
SU(3) symmetry. In Section 5, the dimerized phase with mod2 translation symmetry 
breaking is readily identified by the even-odd asymmetry of the ground state bipartition. 
At the negative biquadratic point, breaking of external symmetries, i.e. symmetry 
transformations involving the spatial coordinates such as translation, together with 
internal symmetry enhancement is shown to yield a rich entanglement structure. In 
the symmetry-protected topological ffaldane phase discussed in Section 6, we find 
a correspondence between entanglement spectra and the low-energy spectrum of an 
open boundary Hamiltonian. We further employ the single mode approximation to 
obtain perturbative ground states around the AfHeck-Kennedy-Lieb-Tasaki (AKLT) 
point, and show that it correctly yields the lowest-lying entanglement deformations 
away from AKLT. In Section 7, we conclude that entanglement spectra for finite size 
realizations constitute a highly efficient tool to understand the underlying structure of 
one-dimensional quantum many-body systems. 

2. Model and ground state bipartition 

We consider the bilinear biquadratic spin 1 chain given by 

H = ^ cos 0 S n ■ S n+ i + sin 6 (S n • S n+ i) 2 , (1) 

n 

where S n = are SO(3) spin operators [S'", Sf,] = ihe a gjS^S nn > and the 

angular parameter 0 < 9 < 2n specifies the sign and relative strength of the bilinear 
and biquadratic interactions. The phase diagram has been thoroughly investigated by 
numerical studies (Fig. 1). A ferromagnetic phase for 7r/2 < 6 < 57t/4 exhibits spin 
rotation symmetry breaking and is the only phase of (1) in which the ground state does 
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Figure 1. (a) Schematic phase diagram of the bilinear biquadratic model (1). The four 
points highlighted by red circles are of central interest in the following: (1) Heisenberg 
ferromagnet, (2) Uimin-Lai-Sutherland point, (3) the negative biquadratic model, and 
(4) the AKLT point, (b) Spatial bipartition of an open boundary 12 site chain resulting 
in one single cut. (c) Spatial bipartition of a periodic 12 site chain resulting in two 
cuts. 

not reside in the singlet sector. The phase for 7t/4 < 6 < 7t/2 has been investigated 
numerically [19] and was found to be gapless with soft modes at k — 0, 2tt/3 and 47t/3, 
with the dominant quadrupolar correlations [20, 21], At 9 = 7t/4, the quadrupolar 
phase terminates in the integrable Uimin-Lai-Sutherland point [22, 23, 24], which will be 
discussed further below. Throughout —7t/4 < 6 < 7t/4, the ground state is a symmetry- 
protected topological phase identified by Haldane [25]. Anchored by the exact solution 
for a gapped valence bond solid at the AKLT point tan# = \ [26], the persistence of 
the spin gap throughout the Haldane phase was investigated and confirmed by early 
large scale numerics [27, 28]. Finally, the dimerized phase for 57t/4 < 6 < 77t/4 exhibits 
translation symmetry breaking as well as an associated dimerization gap, and is framed 
by the multicritical SU(2) 2 Takhtajan-Babujian point [29, 30, 31] to the Haldane phase 
and by the negative Uimin-Lai-Sutherland point to the ferromagnetic phase [32], where 
the latter transition point is challenging to resolve numerically because of the rapidly 
vanishing dimerization gap [33]. 

Consider a finite size realization of (1) for a given 6* for N sites on open (OBC) or 
periodic (PBC) boundary conditions (Fig. 1). Because of finite size splitting, there will 
always be a unique ground state | 'Bq* ) located in a Hilbert space H = lxlxl... xl = 
<8)^1, where 1 represents the local spin-1 degree of freedom at a single site. According 
to the Schmidt decomposition with respect to a spatial bipartition, we split the Hilbert 
space into two regions A and B according to H = hi a x He, with the number of sites 
included in the regions given by Na+Nb = N. From there, we define the reduced density 
matrix p A =Tr bP, where p = I'Lq*) (^o*|; an d T tb denotes a partial trace over all Nb 
sites. Following [2], this is associated with a Boltzmann factor of unit temperature 
p A = exp(— Ha), where Ha denotes the entanglement Hamiltonian of region A. 
Depending on OBC or PBC, Ha mimics one or two cuts between region A and B (Fig. 1). 
The entanglement spectral levels will be denoted by f = spec[i74] = — log(spec[p j4 ]). 
(Note that aside from the real space partition introduced above, other cuts such as the 
momentum cut [34, 35, 36] or the rung cut [37, 38] have been advocated in the context 
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of low-dimensional spin systems.) 

In the case of an isotropic spin chain as in (1), the block diagonal structure of p A 
(and equivalently Ha) is fully determined by the spin representation of |Tq*). If it is a 
singlet, i.e. S a | v b(f ) = 0, it yields [S^, p A ] = 0 , where S A = Yhi&A a — x i Vi z -> which 
gives the SO(3) multiplet structure of [S^,p A ] = 0 characterized by the eigenvalues 
Sa(Sa + 1)- If I'J'o*) i s a s Pi n multiplet such as in the ferromagnetic phase, it only holds 
[S A ,p A ] = 0 which gives a U(l) symmetry rather than a full SO(3) multiplet structure. 
In what follows, N A and Sa characterize the subblocks of Ha- These quantum numbers 
allow to resolve manifestations of symmetry breaking, enhanced internal symmetry, and 
symmetry-protected character from the sole knowledge of finite size ground states of (1). 

3. Ferromagnetic phase: spin rotation symmetry breaking 

The ferromagnet extends to a wide region in (1), hosting at its boundaries at 9 = 7t/ 2 and 
9 = 57t/ 4 two SU(3) symmetric ferromagnets. In many respects, the ferromagnet is the 
most trivial phase from the viewpoint of entanglement. In the thermodynamic limit were 
a magnetization axis is spontaneously chosen, the ground state is a single site product 
state and does not contain any entanglement. For finite size, however, the continuous 
spin symmetry of (1) cannot be broken, whereas the ground state I^fm) picks the 
largest available SO(3) representation. As such, |T F m) is uniquely defined independent 
of the specific choice of Hamiltonian in the ferromagnetic phase. As elaborated on in 
Section 2, |T F m) retains a U(l) symmetry in the entanglement spectrum p A such that 
we can decompose it into according subblocks. We restrict the following discussion to 
S z > 0 (S A > 0) as the U(l) symmetry guarantees that subblocks with —S z (—S A ) are 
identical to those with S z ( S A ). 

We begin with the ferromagnetic spin 1 ground state wave function. In the spin-p 
case, the wavefunction has the most simple form, as all basis states acquire identical 
weight. For spin 1, we start by the subsector with maximum magnetization S z = N, 
= ITT - - - T)• From the knowledge of the multiplet structure, it must hold 
S 2 0n m |*Im =S > = S (S + 1 )(^“) m |^lM S > Vm G [0,...,2S]. This allows us 
to compute the relative weights of the different wave function partitions in S z basis. 
With respect to I^fm""^)) the spin 1 ferromagnetic wave function configurations acquire 
different relative weights due to different multiplicities of spin flip operators acting on 
the same site. For example, given S z = N — 2, the basis is constructed by applying 
S~ Sf [ I f ... t), and configurations with i ^ j acquire twice the weight of i — j due to 
bosonic normalization. 

Consider the bipartition of |Tfm) in the S z = 0 sector with N A = N/2, as illustrated 
in Fig. 1. In the following, we assume the Schmidt basis of region A and B to be 
ordered according to decreasing relative weight. From the Schmidt decomposition, the 
(in general rectangular) entanglement matrix Ea can be defined as p A = E a Ea , which 
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Figure 2. (Left) Analytical entanglement spectrum of the spin 1 ferromagnet as 
a function of total system size N for PBC. Each S\ sector contains one finite 
entanglement level. Different colors correspond to different S A sectors as indicated 
in the legend. (Right) Example: ES vs. S A for N = 12. 


for S Z A = 0, aside from global normalization, takes the form 


E a = 


( 1 1/2 1/2 
1/2 1/4 1/4 
1/2 1/4 1/4 

V ; 


1/4 1/4 
1/8 1/8 
1/8 1/8 


(l/2) m 

(l/2) m+1 

(l/2) m+1 


\ 


( 2 ) 




One finds that the matrix E A consists of several blocks with identical individual entries, 
the size of which being dictated by the number of basis states with the same relative 
weight. Determining the size of these blocks is given by the dimension of the S z (S A ) 
subspaces of the Hilbert space for a given chain length N (Na), 

l(N-S*)/ 2J 


TV 

S z + 2 j 


S z + 2 j 
J 


( 3 ) 


dim % (TV, S Z ) = J2 
j =o 

where denotes the floor function. By construction, all columns (and rows) are 
linearly dependent, i.e., Ea has rank 1 and the characteristic polynomial reads P[/i] = 
H Na ~ 1 (// — g ) where p = g is the only non-zero eigenvalue given by 


L(AA-Si)/2J 
9 Wa,s* a ] = 

3=0 


N a 
+ 2 j 


^ S A + 2 -(2 j+S z A )_ 


( 4 ) 


Restoring the overall normalization of leads to 9[ Na s-] = 9[n a s z ]/V^ N ’ S 'totl 5 

where the norm is given by 


L(7V-5f ot )/2j 
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Figure 3. (N,Na) = (12,6) entanglement spectrum (PBC) of the Uimin-Lai- 
Sutherland point decomposed into SO(3) representations labeled by Sa- The enhanced 
SU(3) symmetry becomes apparent from the eigenvalue degeneracies of different SO(3) 
representations, such as a triplet and quintuplet forming an SU(3) octet as the first 
entanglement energy above the lowest singlet level. 



The non-zero eigenvalue of the reduced density matrix is then given by 


g [n a ,s z a 


, where 


the entanglement spectrum £ = — log \g^ Nji S z j \ is visualized for large system sizes in 


Fig. 2. In the right panel of Fig. 2, the same ES for N = 12 is plotted as a function of 
Sa- The parabolic form is a 
for the different S Z A sectors. 


Sa- The parabolic form is a direct consequence of the binomial Hilbert space dimension 


4. Quadrupolar phase: enhanced SU(3) symmetry at 6 — 7t/4 

Following numerical work [19], it was understood that the quadrupolar phase could be 
suitably characterized by the 6 = 7t/4 point. This is the Uimin-Lai-Sutherland (ULS) 
model, an integrable SU(3)i model with an exact solution of the Bethe Ansatz type 
[22, 23, 24], At this point, the interaction between sites is, up to a constant, proportional 
to — P 1 (n, n + 1), where P x is the projector onto total spin J = 1 for the neighboring 
sites. The J = 1 triplet is thus the ground state of the link, and the J = 0 singlet and 
J = 2 quintuplet combine to form a degenerate sextuplet. The Hamiltonian can then be 
mapped onto that of the antiferromagnetic SU(3) Heisenberg chain in the fundamental 
representation, since 3®3 = 3ffi6. As known from the SU(iV) formulation of the 
Lieb-Schultz Mattis theorem [39], and consistent with the generalized SXJ(N) Haldane- 
gap criterion [40], this model is gapless. Despite its solvability, from the viewpoint of 
numerical entanglement measures, it is still generally involved to extract the critical 
theory content for such models from finite size scaling [41]. 

The enhanced internal SU(3) symmetry, however, can be readily derived (see 
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Figure 4. Young tableaux decomposition of all SU(3) representations contained 
in the Na = 6 Schmidt basis. The representations that do not comply with the 
SU(3) bipartition sum rule for an SU(3) singlet ground state are highlighted in orange, 
corresponding to > 5 symmetrized blocks. 


also Appendix A) and, because of its locality, manifests itself for any finite size 
realization. Along with the SU(3) singlet property of the ground state, this yields an 
SU(3) multiplet structure for p A , where the representations of SU(2) arrange accordingly 
in the entanglement spectrum (Fig. 3). (Alternatively, one can directly adapt SU(3) 
Casimir operators J 3 and J 8 to label the individual blocks of p A .) Aside from the exact 
SU(3) symmetry specific to the ULS point, the qualitative entanglement features carry 
over to the whole quadrupolar phase. 

A notable aspect of p A for finite size realizations of the ULS point is the observation 
of an extensive number of zeroes which are due to the enhanced SU(3) symmetry. To 
illustrate the latter for N = 12, note first that because of its singlet property, any 
ground state configuration should contain 4 green, blue, and red particles according 
to the notation in Fig. Al, which directly connects to 4 sites located in the state 
Sf = —1, 0, and 1, respectively. Considering now the Schmidt basis of region A 
(Fig. lb), it trivially follows that any basis configuration will contribute a zero eigenvalue 
which would necessarily yield a total state configuration violating the SU(3) sum rule 
mentioned before. For IV = 12 and Na = 6, this corresponds to all configurations with 
5 or 6 particles of one color in the Schmidt basis. Since all configurations form SU(3) 
multiplets, such configurations are interpreted as the seed states for the multiplets they 
belong to which are all disallowed, i.e. the total number of zeroes is given by the total 
dimensionality of the associated SU(3) representations. How many such representations, 
and hence how many zeroes of p A , there are can be elegantly expressed in terms of SU(3) 
Young tableaux. For N — 3 * n, n G N, it corresponds to counting all representations 
with a number of symmetrized blocks that exceed n (Fig. 4). In our case, we have to 
consider 3 06 = 5 • 1 © 16 ■ 8 © 10 • 10 © 5 • 10 © 9 • 27 © 28 © 5 • 35, and the total 
number of zeroes hence gives 5 x 35 + 28 = 203 in the ground state bipartition p A due to 
SU(3) symmetry. Away from the ULS point, the zero modes of p A for ULS become very 
small eigenvalues of p A and hence yield entanglement levels at high £. It gives a generic 
explanation for the large amount of high entanglement energy levels in the quadrupolar 
phase. The argument of bipartition-induced zero mode representations carries over to 
analogous scenarios for SU(iV) entanglement spectra. 
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Figure 5. N = 12 entanglement cut for Na = 5 and Na = 6 for the dimerized SU(3) 
point at 0 = 3ir/2 for PBC. Aside from degeneracies between different Sa sectors 
dictated by enhanced SU(3) symmetry, the low-lying entanglement spectrum of the 
gapped state switches between a singlet and triplet, which is reminiscent of translation 
symmetry breaking in the dimerized phase. 


N A =Q 


3 Sa 4 



5. Dimerized phase: Translation symmetry breaking 

5.1. SU(3) dimerized point at 6 = 37t/2 

As alluded to in Appendix A, the model (1) for 6 = 3n/2 exhibits enhanced internal 
SU(3) symmetry similar to the ULS point discussed in Section 4. Here, however, the 
Hilbert space relates to an alternating 3 and 3 representation. This has fundamental 
consequences on the very nature of the state. While it is still amenable to analytic 
solution, the state is not gapless but exhibits a gap along with translation symmetry 
breaking [42, 43], where the gap could be quantified by transfer matrix Bethe Ansatz [44] 
and allows to estimate a spin-spin correlation length of 21 sites. It is revealing how the 
joint appearance of enhanced SU(3) and translational symmetry breaking manifests 
itself in the entanglement spectrum for small system sizes. To begin with, because of 
the 3x3 Hilbert space structure, there is only an SU(3) sum rule for neighboring sites: in 
terms of SO(3) spin configurations for an even number of lattice sites, it implies that one 
finds non-zero weights in the ground state wave function only for basis states composed 
of two-site sequences [S'f,S'f +1 ] = [1,-1], [—1,1], or [0,0]. As such, these configuration 
constraints only act locally on consecutive sites, and do not trigger extensive zero modes 
as observed in Section 4. 

5.2. Entanglement gap vs. energy gap 

The onset of translational symmetry breaking can be detected e.g. by analyzing 
(. N,N a ) = (12,5) and (12,6) (Fig. 5). For gapped spin chains, the correspondence 
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Figure 6. Upper panel: spectral flow of entanglement spectra in the dimerized 
phase obtained within DMRG. Lower panel: energy gap (green) in units of (1) and 
entanglement gap (red) in units of Boltzmann temperature for N = 48 and N = 120 
(OBC). The entanglement gap maximum is shifted against the energy gap maximum. 


between the entanglement spectrum and the open boundary Hamiltonian spectrum can 
be developed (see e.g. Ref. [45]). This is because the low-energy modes contributing to 
the entanglement Hamiltonian, as located at the Schmidt boundary , correspond to the 
low-energy modes located at the physical boundary. For entanglement spectra in the 
dimerized phase, the lowest lying singlet state seen for Na = 0 mod 2 alternates with the 
lowest lying triplet state for Na = 1 mod 2. Note that this is not trivially dictated by 
the Hilbert space structure, as any even or odd tensor product of SO(3) representations 
contains singlets. Upon inspection, both low-energy spectra for (TV, Na) = (12, 5) and 
(12,6), in terms of spectral structure and eigenvectors, correspond to an OBC TV = 5 
and TV = 6 realization of (1). The even chain yields a lowest-lying dimerized singlet, 
while the lowest-lying odd chain state breaks a dimer singlet into a triplet state. 

This spectral structure carries over to the entanglement spectrum of the whole 
dimer phase. For the Na = 6 cut, an entanglement gap between the singlet state 
and the second lowest quintuplet state ( 6 > 37t/2) or triplet state (6 < 37t/2) can 
be defined (c/. Fig. 6), with a corresponding level crossing cusp at 6 = 3tt/2 [46] 
(see also Appendix B). Interestingly, the energy gap shows a similar, but not identical 
behavior. As can be also supported by larger scale density matrix renormalization group 
(DMRG) calculations, the energy gap maximum does not match the entanglement gap 
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Figure 7. Entanglement spectrum for (N,Na) = (12,6) and (12,5) in the Haldane 
phase at 0 = 7r/20 for PBC. The low energy spectrum is given by a triplet and a 
singlet separated from higher levels. The spectrum is similar to the corresponding 
Hamiltonian spectrum for open boundaries. 


maximum (Fig. 6). It stresses once again that while certain universality features are 
shared between energy spectra and entanglement spectra in gapped phases, this does 
not extend to non-nniversal features such as gap maxima. 

6. Symmetry-protected topological Haldane phase 

The Haldane phase exhibits an energy gap without breaking of translation symmetry. 
The entanglement analysis shows no breaking of any other symmetry of (1). The 
correspondence between the entanglement spectrum and the OBC Hamiltonian 
spectrum can be studied already for small finite system size. Fig. 7 shows the 
entanglement spectrum for even and odd N&. The general structure of the entanglement 
spectrum stays unchanged, in clear contrast to the analogous analysis for the dimerized 
phase in Fig. 5. A triplet and a singlet level appear to form a separated low-energy 
set from the rest of the energy spectrum, where the singlet (triplet) state is the lowest 
entanglement energy state for even (odd) Na- This is identically found for the OBC 
Hamiltonian spectrum. 

6.1. Entanglement gap vs. energy gap 

The gap of the Haldane phase, and its entanglement Hamiltonian correspondence to 
OBC spectra derived from there, is identified as a Berry phase effect in integer spin 
chains [25] and characterized by the AKLT point where the gap can be calculated 
analytically [26]. The AKLT model has become the paradigmatic symmetry-protected 
topological phase in one dimension [11], where the OBC low-energy behavior is 
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Figure 8. Upper panel: Entanglement spectral flow of the Haldane phase obtained 
within DMRG for IV = 48 (PBC). Lower panel: Energy gap vs. entanglement gap 
for OBC. The proper investigation of the energy gap for large system sizes yields 
gap closure at ±7r/4 and a maximum energy-gap position shifted against the AKLT 
point [28]. 


characterized by one dangling spin-| degree of freedom at each boundary. This likewise 
characterizes the low-energy entanglement spectrum, i.e., the dangling spins form a 
triplet and a singlet which are separated from the other levels by an entanglement gap 
or energy gap, respectively (Appendix B). While the energy gap always stays finite, the 
entanglement gap becomes infinite at the AKLT point, where the ground state can be 
exactly characterized by a matrix product state with bond dimension V = 2. 

Fig. 8 depicts the entanglement spectral flow through the Haldane phase along with 
a comparison of energy gap vs. entanglement gap. From the ALKT point, branches of 
levels come down in entanglement energy and eventually close the entanglement gap 
towards the ULS point at 7t/4 and the Takhtajan-Babujian point at —vr/4, respectively. 
The energy gap similarly closes at ±7t/4 as seen in finite size scaling, but exhibits a 
gap maximum shifted against the entanglement gap maximum at the AKLT point. 
(The observation from small finite size is confirmed by large scale calculations, see e.g. 
Ref. [28].) This is similar to the dimerized phase in Fig. 6. 
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Our starting point is the AKLT Hamiltonian plus the Heisenberg term as a small 
perturbation, 

# = tfAKLT + ^S;-S m . (6) 

i 

Using the spectrum of single mode approximation (SMA) excitations, it is possible to 
derive an approximate perturbation theory in the difference H— Haklt- (Note that SMA 
was employed previously to elucidate the relation between valence bond states (VBS) 
and the Laughlin state [47].) The idea goes back to unpublished work of Haldane, and 
was later elucidated in Ref. [48]. One proceeds by idealizing the triplet SMA levels as 
bosonic excitations, and asserting the approximate operator correspondence 

Si « VW) (*U + b- t J, (7) 

where s(k) = (T 0 |5'fc<S'“ i .|T 0 ) (no sum on a) is the static structure factor. The bosonic 
form of the perturbed Hamiltonian in (6) is then 

H hos = J2 + £s ( k ) cos k{b[ Q + b_ Ka ){b ] _ Ka + b ka ) , (8) 

k,a 

where the structure factor and SMA dispersion were computed in Ref. [47] and found 
to be s(k) = 2(1 — cos/c)/(5 + 3 cos k) and co(k) = A(5 + 3 cos k). Solving H hos by a 
Bogoliubov transformation, one obtains the renormalized structure factor 

™= LvTLw Y ™' (9) 

where A (k) = 2 s(k) cos k. The corresponding spin wavefunction is then given by 


T 0 \ = exp - ^ 


s{k) s{k) J SkS ~ k I ^ ’ 


( 10 ) 


where |T 0 ) denotes the AKLT ground state. Expanding to leading power in e, we have 


T n ) = 


v n,7i' J 


^o) , 


( 11 ) 


where the kernel K nn > is the Fourier transform of 2A(k)/oj(k)s(k) = 4 cos k/oj(k), and 
is given by the expression (j := n — n') 

2tt 

lh Inns lc ... .“tfi / 1 \ 3 

( 12 ) 


K ^ J 27t u(k) 


dk 4 cos k iki 36 / ly' 


The expression for 
wherein 


Tq) has a clear correspondence with first order perturbation theory, 


( 13 ) 


f'o) = wm 2.y cos (k)sisi t n„>, 


if we approximate | n) ~ s(k) 1 (1 — |T 0 ) (T 0 |) S k Sf k |T 0 ), and E n ~ c u(k). 
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Defining the operator-valued matrix, 

/ 1 

0 

M(j) = ' 




S] 

0 


0 
0 
Vo 


0 

0 

0 


Sj 

0 

0 


o \ 


0 

0 


Sf 

S] 

Sf 


(14) 


and the vectors 


(L|=( 


-3e 0 0 0 1 


|R) = 


/ 0 \ 
0 
0 
0 

V 1 / 


(15) 


we have the matrix product operator (MPO) expression, 
(L| M( 1) • • • M(N) |R) = 1 + K ( n 

n<n' 


n 


QCt QQ! 


(16) 


If the AKLT state is written in matrix product form as 

|'ho) = {C\A mi ■ ■ -A mN \1Z) \mi, ..., m N ) , (17) 

where (£| and (771 are vectors which passivate the end sites n — 1 and n = N, rendering 
hen the SMA state is a MPS with 

\V) = (L\A m '---A m » \R) | m u ...,m N ) , (18) 


them S 

= then 


m 

where 



Am 

A i,i' 

and i = 

(a, /i) is a 




(19) 


and m from 1 to 3. Thus, i runs from 1 to 5 • 2 = 10. Similarly, | L) = |L) <g> |£), etc. 
The fact that A is of rank ten means that there must be ten entanglement levels. 

More frequent than a double cut of a finite size bipartition, the single cut associated 
with a semi-infinite partition has been employed to identify the symmetry-protected 
topological character of the Haldane phase [11]. The latter is simulated for small finite 
size by passivating the dangling spins in an OBC geometry (Fig. lb) and choosing the 
S* ot = ±1 sector. Combining the valence-bond picture (see Fig. 9) and the framework 
of SMA, we will demonstrate that the ES in the Haldane phase must be two-fold (or, 
in general, even-numbered) degenerate. 

Suppose we partition the chain such that sites j G {1,2, ...,77a} are in subsystem 
A and the remaining sites are in subsystem B. It follows that we may write 

I^sma) = l^o) + I^a) + |^b) + I^ab) , (20) 

where |4 >a) includes contributions for n < n' < N\, |Tb) includes contributions for 
TVa < n < n', and (Tab) includes contributions for n < N& < n' (Eq. 11). 
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(a) 




12 n n! N—1N 


N -1 iV 


Figure 9. (a) S = 1 AKLT state; each lattice site (pink ellipses) is given by 

the symmetrization of two virtual spin-^ degrees of freedom. One of them is 
antisymmetrized in the singlet bond to the left while the other is antisymmetrized 
in a singlet bond to the right, (b) Visualization of the state | nn') within the SMA 
framework for the AKLT chain with passivated S = 1 ends. 


We may then write the density matrix p = |\&) (4/| as p = p 0 + p A + p B + p AB , 
where p A = |'P 0 ) (V A \ + |^a) (^o|, Pb = l^o) (^b| + |^b) (^o|, and p AB = |tf 0 ) (^ab| + 
|4>ab) (4/o|. We now perform the partial trace over the B subsystem. For p 0 , p A , and p B , 
this results in one link being severed by the entanglement cut. There are then exactly two 
entanglement eigenstates in each case (he., a doublet). For p AB , three links are severed 
(see Fig. 9), resulting in two doublets and one quadruplet entanglement levels. Assuming 
that the doublets are linearly independent, we find three linearly independent doublets 
and one quadruplet, resulting in ten entanglement levels, as previously deduced. The 
upshot of this analysis is that one predicts the appearance of doublets and quadruplets, 
but nothing else. Here we recover the two-fold degeneracy (or, more generally, even- 
number degeneracy) of the ES for the Haldane phase [11]. If we allow for higher processes 
and go beyond the linearized SMA version considered here, we will generate higher-order 
links as compared to | nn'). Still only configurations are generated where the real space 
cut intersects an odd number of links. Therefore the ES will remain even-numbered 
degenerate. 

The double degeneracy of the ES for a single cut in the Haldane phase clearly 
distinguishes the Haldane phase from the dimer phase. While the SMA picture can 
explain the two-fold degeneracy in the ES, one might naively expect that breaking of 
SU(2) symmetry would lift these degeneracies. This is, however, not true. As shown in 
Ref. [11], the Haldane phase is protected by time-reversal, bond-inversion, and dihedral 
symmetries. One needs to break all these symmetries in order to loose the degeneracy 
of the entanglement spectrum. In the meantime it is well-established that the Haldane 
phase is an SPT phase [49] . The double degeneracy in the entanglement spectrum is a 
hallmark for such a SPT phase when only a single cut is considered, as just explained. 
When PBCs are imposed and two cuts are present, multiplets stemming from the two 
cuts decompose into multiplets such as singlet, triplet etc. The simplest example is the 
pure AKLT state where a doublet on both cuts is present resulting in one singlet and 
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Figure 10. Comparison of the SMA-corrected ground state entanglement spectrum 
(leading order in e, blue) (Eq. 11) and the exact ground state entanglement spectrum 
for a single cut in a IV = 12 site chain of the Hamiltonian in Eq. 6 (red). All 
entanglement levels are even-fold degenerate in both cases. Since we mimic the single 
cut by imposing OBC and S* ot = 1, exact SU(2) symmetry in the ES is lost, but is 
still approximately present. 


one triplet level. There is no non-trivial even/odd degeneracy constraint anymore. 

6-4- Comparison of SMA predictions with numerical results 

We test the SMA method quantitatively by comparing the entanglement spectrum of 
the SMA-modified ground state wavefunction with that of the exact ground state. We 
first consider an open chain with iV = 12 sites and make a single entanglement cut which 
divides the system into two half-chains. The end spins are passivated by fixing S* ot = 1. 
In Fig. 10 we show both SMA (blue) and exact ES (red) for e = 0.001. The lowest- 
lying doublet, the entanglement ground state, is followed by a quadruplet and another 
doublet. The SMA state entanglement level multiplet structure is as predicted above 
(i.e. ten levels in total), but deviates noticeably from the exact result for entanglement 
levels above the lowest doublet. One reason is that the SMA-corrected ground state is 
of the form exp e 3 Qj'j |T 0 ), where {Qj} are operators which scale extensively 

with system size. Hence an expansion of the exponential is really an expansion in 
powers of Ne. We can do a little better by expanding the SMA-corrected ground state 
wavefunction to second order, writing 

®o) = {l - |x (1) + j X (2) + J (M m ) 2 + d(e 3 )| |f„) , 


( 21 ) 
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Figure 11. Evolution of entanglement spectra for an N = 12 site chain of Eq. (6) as a 
function of e. A single cut is used with S* ot = 1 enforced to passivate the two end spins. 
Comparisons of the exact ground state (red) with first (blue) and second (green) order 
SMA-corrected ground states are shown. All states exhibit the characteristic AKLT 
doublet whose £ level does not change with e within machine precision. All doublets 
are depicted by full circles and the quartets by full squares. The lowest entanglement 
levels above the AKLT doublet are hardly affected by moving from 1st to 2nd order 
SMA. The quantitative deviation of the SMA levels from the exact levels is small in 
terms of absolute entanglement weight, given the high values of £. 


where XW = Er A) S ? S ~* 

Eq. 12 above, and 


E 


7 >'( r ) Ca oa 
n,n' ^n' 


with hf ] 



864 27 

~25~ dj ’° + EO 


(~l Y- (527 — 180j + 400j 2 ) . 


K{j) in 
( 22 ) 


In Fig. 11 we compare the low-lying entanglement levels related to Eq. 6 for the exact 
ground state as well as the first and second order SMA results, as a function of e. We 
End that the ES is strongly affected for levels above the lowest eight states arranged 
in the doublet, quadruplet, doublet order. Due to the second order contribution, 
the sequence of entanglement levels (d=doublet and q=quartet) changes from d-q-d- 
d for the first order SMA to d-q-d-q, which matches the exact sequence. Still, the 
second order contribution modifies the entanglement profile only marginally. The fact 
that there is still only a rough quantitative correspondence in £ between the higher 
lying entanglement levels for the SMA-corrected and exact ground state indicates that 
the source of the disagreement lies in the approximate nature of the Bogoliubov-SMA 
approach itself, rather than the passivation of the end spins by projecting onto Sf ot = 1, 
or the finite order expansion in Ne. 
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We have employed finite size ground state bipartitions to analyze the entanglement 
structure of isotropic spin-1 chains. The isotropic bilinear-biquadratic model allowed 
us to interpret many hallmark features of quantum spin chains from the viewpoint 
of entanglement. We investigated spontaneous breaking of translation symmetry in 
the dimerized phase and spin rotation symmetry in the ferromagnetic domain. Both 
phenomena can be unambiguously identified from entanglement just by analyzing 
a single small system size. Furthermore, we have shown how enhanced internal 
SU(3) symmetry affects not only the degeneracy structure of entanglement spectra 
similar to Hamiltonian spectra, but also the pronounced extensive amount of zero 
weight entanglement levels as seen for the Uimin-Lai-Sutherland point. This spectral 
distribution feature of entanglement approximately persists for the whole quadrupolar 
phase. Finally, the SPT character of the Haldane phase is precisely resolved by 
entanglement spectra. We have elucidated the correspondence between a boundary 
Hamiltonian spectrum and the entanglement spectrum of the associated ground state 
for periodic boundary conditions. In particular, we have identified the notion of 
the entanglement gap to hold in the thermodynamic limit. For the single mode 
approximation of the perturbed AKLT state and the operator product expansion derived 
form there, entanglement spectra allowed us to obtain a complementary view on the 
accuracy of this approach. 

Our study strongly supports the view that the analysis of entanglement spectra 
of one-dimensional quantum systems provides insight which leverages the analysis of 
relatively small finite size systems, and will thereby constitute a preferable choice when 
large scale calculations are either not feasible or unavailable. 
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Appendix A. Derivation of enhanced SU(3) symmetry models 

For completeness, we elaborate on the points of (1) where the Spin 1 SO(3) model 
exhibits enhanced internal SU(3) symmetry. We follow the notation of [50], and prove 
below that both Hamiltonians 6 = 7t/4 and 6 = 3n/2 of (1) possess an enlarged SU(3) 
symmetry. Since both models exhibit SU(3) singlet ground states, the same argument of 
for SU(2) in Section 2 applies and as such, the internal symmetry operator generating the 
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block diagonal form of the entanglement spectrum is given by [T a , p A ] = 0, T a = JA Tf, 
where T a are the generators in sl(3). This yields SU(3) multiplets in the entanglement 
spectrum. 

Appendix A.l. SU(3) Heisenberg model 

To prove the SU(3) invariance at the Hamiltonian level, we develop an a posteriori 
perspective and consider the SU (3) Heisenberg Hamilton operator 

N 

(A.l) 


H = Y. T - T w 


with implicit summation over a = 1,... 8. We label the representation states of the 
fundamental representation rep. 3 by the colors blue (b), red (r), and green (g) 
(quarks). We will also treat the non-equivalent representation 3 later, i.e., where the 
states possess the complementary colours yellow (y), cyan (c), and magenta (m) (anti¬ 
quarks). Labelling the Gell-Mann representation by A“ in order not to confuse it with 
the generators T, the action on the fundamental representation 3 is given by 
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(A.2) 


0 0 0 \ 

0 0 -i 

0 i 0 

Their normalization is chosen to be 
tr (A a A 6 ) = 2 5 ab . 

The Gell-Mann matrices form an orthogonal basis of sl 3 , the Lie algebra of SU(3), and 
satisfy the commutation relations 

[A°, \ b ] = 2if abc X c . (A.3) 

The structure constants f abc are totally antisymmetric and obey Jacobi’s identity 

j^abc jcde _|_ jtbdcjtcae _|_ jtdac jtcbe q 

All non-vanishing structure constants are obtained by permutations of the indices from 
/ 123 = 1 


j 147 __Z^ 40 _— 


f 246 


257 


f345 


j 156 ^36Y ' 


367 


f458 _ j 678 . = ^3 


(A.5) 
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Figure Al. Weight diagrams of the three-dimensional representations of SU(3). T 3 
and T 8 are the diagonal generators. 

The Gell-Mann matrices also close with respect to the anticommutator: 

4 

{A a , ^b } T d a \) C \ C i 

with the fully symmetric structure constants 

d 118 = d 228 =d 338 = -d 888 = 4= 

V3 

d 448 = d 558 = d 668 = rf 778 = _J_ 

V3 

d 146 = d 157 = -d 2A7 = d 256 = h 344 = d 355 = -d 3m = -d 377 = 1 

Appendix A.2. SU(3) symmetry at 6 = 7 t /4 

We now specify the Hilbert space the Hamiltonian (A.l) acts on. Considering 
3 x 3 x 3 x 3 x 3..., we can thus write the Hamiltonian more specifically in terms 
of the As: 

N 

ff = E A ? A ?+1- ( A ' 8 > 

i 

Next, note that the fundamental representation of the Lie algebra o(3) relates to the 
SU(2) spin 1 operators, because the SU{ 2) is locally isomorphic to 0(3) and the 
representations have the same dimensionality. As such, we can readily identify the 
0(3) type operators of the Gell Mann matrices. For s/, Z = 1,2,3 and spin 1 operators, 
this yields 

si = A7, S2 = — A5, s 3 = A2, (A. 9 ) 

enabling us to write the spin 1 bilinear terms in terms of As. The same can be done for 
the biquadratic terms, exploiting anti-commutation and commutation relations of the 
Gell-Mann matrices A [50]. Accordingly, the Hamiltonian part Sj • Sj+i is readily 
identified via A 2 = S3, A 5 = — £ 2 , A 7 = Si and yields 

N N 

^s,-s w = Yi E A ? A “«- 

i i a= 2 , 5,7 


(A.6) 


(A.7) 


(A.10) 
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In turn, the biquadratic Hamiltonian part Hb q = (S,; ■ Sj + i) 2 is rewritten as 


= (A^ +1 ) 2 + (\ 5 AU 2 + (AJA[ +1 ) 2 

(A. 11) 

+ A 2 A'/A 2 ; | A/., + [2 -H- 5] 

(A. 12) 

+ A 2 AjA 2 +1 Aj +1 + [2 o 7] 

(A. 13) 

+ A; 5 A 7 A 3 +1 A 7 +1 + [5^7] 

(A. 14) 


Whenever spatial indices are suppressed in the following, we assume that multiple powers 
of A act on the same site. To calculate Eq. A. 11, we make use of A.6 and find (A 2 ) 2 = 
§(§ + ^A 8 ), (A 5 ) 2 = |(| - T=a 8 + A 3 ), (A 7 ) 2 = |(| - ^A 8 - A 3 ). Factoring everything 
out, Eq. A.11 becomes |A 8 Af +1 + |A 3 A 3 +1 + |. The latter constant we can discard as a 
global constant factor, the rest are just the bilinears of the Casimirs. The subsequent 
lines A.12, A.13, and A.14, can be calculated analogously, so we will only show 
Eq. A. 12. Trivial algebra yields the (anti-)symmetrized form A 2 Af A 2 +1 Af +1 + [2 5] = 

|[A 2 , A 5 ]j{A 2 , A J }j + i + ^[A 5 , A 2 ]j{Ay A J }, + i +1[A 2 , A 5 ]j[Ay A 5 ]j + i + ^{A 2 , A 5 }*{A 2 , A 5 }i+i = 
|[A 2 , A 5 ]j[A 2 , A 5 ]; + i + |{A 2 , A 5 },{A 2 , A 5 };+i. Now, we can make use of (A. 3) and (A.6) to 
End [A 2 , A 5 ] = i\‘ and {A 2 , A 5 } = A 6 , such that the final result is |(A®A® +1 — A 7 A 7 +1 ). 
We End 

N „ N 


X> ■ s i+ o 2 =. E 


v,a=l,3,4,6,8 


\ a \ a 
\ \+1 


a= 2,5,7 


\a \ a 
\ \+1 


(A. 15) 


Adding (A. 10) and (A. 15) yields the Hamiltonian to be symmetric with respect to the 
Gell-Mann index a, i.e. to SU(3). 


Appendix A. 3. SU(3) symmetry at 9 = 3ir/ 2 

Again starting from A.l, we consider the Hilbert space spanned by 3 x 3 x 3 x 3. This 

has consequences for the SU(3) Hamiltonian, as the imaginary-valued A’s for a = 2, 5, 7 
acting on 3 yield one minus sign due to complex conjugation. Performing the same 
algebra as before, we find the purely biquadratic Hamiltonian (A. 15) to be SU(3) 
invariant if acting on the modified Hilbert space, anticipating a doubling of the unit 
cell along the dimerized ground state at 9 = 37t/2. 


Appendix B. Scaling of the entanglement gap 

As discussed in Sec. 5, the even-odd behavior of the ES with respect to entanglement 
cuts corresponding to even and odd subsystem lengths is reminiscent of translation 
symmetry breaking in the dimerized phase. Different, even though similar even-odd 
effects, however, can also be observed under different circumstances such as for for a 
gapless spin- \ Heisenberg chain. For the even-odd discrepancy to provide substantiated 
information, it is hence important to distinguish a gapped from a gapless phase. It 
turns out that the scaling of the entanglement gap clearly distinguishes both cases 
(we restrict the discussion to half-chain cuts). We define the entanglement gap as the 
gap between the entanglement levels belonging to the entanglement ground state and 
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Figure Bl. (Color online) Scaling of the entanglement gap for gapped and gapless 
models for OBC and PBC obtained within DMRG. While the gapless models exhibit 
non-trivial spectral flow for large system sizes (see text for details), the gapped models 
show saturation of the entanglement gap. 


higher-lying entanglement levels. For instance, for spin 1 Heisenberg chain with PBC 
the entanglement gap is above the lowest singlet and triplet; for the spin-f Heisenberg 
chain, we consider the gap above the single lowest singlet state. Fig. Bl shows examples 
for gapped and gapless models. Plots on the left (right) correspond to OBC (PBC). 
The spin 1 Takhtajan-Babnjian (TB) model (6 = T7T/4), the spin-| Heisenberg chain, 
and the spin-f XXZ chain with Ising anisotropy A = 0.5 are examples for the gapless 
case, the spin 1 Heisenberg chain (6 = 0) and the dimer model with 6 = 3n/2 are 
examples for the gapped case. For the spin-), Heisenberg model and the spin 1 TB 
model, both £i — £o and £2 — £1 are shown. The characteristic scaling behavior vs. 1/N 
agrees with the prediction of Calabrese and Lefevre [51, 52], For gapped models, we 
find that the entanglement gap already saturates for rather small systems sizes (N = 30 
for OBC, N — 60 for PBC), supporting the existence of an entanglement gap in the 
thermodynamic limit. 
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